UV Completions of Magnetic Inelastic Dark Matter and RayDM for the Fermi Line(s) 



Neal Werner 1 ^ and Itay Yavin 2,3 '[j] 

1 Center for Cosmology and Particle Physics, Department of Physics, New York University, New York, NY 10003 
2 Department of Physics & Astronomy, McMaster University 1280 Main St. W. Hamilton, Ontario, Canada, L8S 4L8 
3 Perimeter Institute for Theoretical Physics 31 Caroline St. N, Waterloo, Ontario, Canada N2L 2Y5. 

Models that seek to produce a line at ~130 GeV as possibly present in the Fermi data face 
a number of phenomenological hurdles, not the least of which is achieving the high cross section 
into 77 required. A simple explanation is a fermionic dark matter particle that couples to photons 
through loops of charged messengers. We study the size of the dimension 5 dipole (for a pseudo-Dirac 
state) and dimension 7 Rayleigh operators in such a model, including all higher order corrections 
in 1/M mess . Such corrections tend to enhance the annihilation rates beyond the naive effective 
operators. We find that while freezeout is generally dominated by the dipole, the present day 
gamma-ray signatures are dominated by the Rayleigh operator, except at the most strongly coupled 
points, motivating a hybrid approach. With this, the Magnetic inelastic Dark Matter scenario 
provides a successful explanation of the lines at only moderately strong coupling. We also consider 
the pure Majorana WIMP, where both freezeout and the Fermi lines can be explained, but only at 
very strong coupling with light (~ 200 — 300 GeV) messengers. In both cases there is no conflict 
with non-observation of continuum photons. 



I. INTRODUCTION 

One of the most striking signals for indirect detection 
of dark matter is a monoenergetic gamma ray line. Such 
a feature, in different locations, with the same energy 
has no known astrophysical background. While by its 
very nature dark matter must only couple very weakly 
to photons if at all, even a small signal could be very 
convincing. In light of this, the recent claims of a line 
from the galactic center at approximately 130 GeV [IHl] 
(and possibly a second near 111 GeV[l]) are extremely 
exciting. Subsequent and statistically weaker claims to 
have seen a consistent signal in galaxy clusters [5] and 
unassociated point sources [5] (although see 0(5]), while 
individually perhaps unpersuasive, make the possibility 
even more intriguing. 

There are possible concerns, of course. The location 
of the peak in the GC seems not precisely central [4] (al- 
though see [S]), contrary to naive expectations. However, 
some numerical simulations that accounts for both dark 
matter as well as baryons have found possible preliminary 
evidence for such a transient feature [10]. No assoicated 
continuum emission has been detected from the galactic 
center [TlTU3| . Diffuse limits provide tensions, but no 
clear exclusions [2]. A strange peak at the same energy 
has been found in a subset of the Earth limb photons 
|15j , raising some concern that a bizarre instrumental ef- 
fect could be giving rise to this. At the same time, no 
compelling explanation has been offered that would yield 
a sharp peak in the GC and not in e.g., the disk as well. 

Regardless of these issues, it is worth understanding 
what models could explain these data. A number have 
been put forward (e.g., [T6"rl2"T] ) . An effective theory 
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approach was taken in 20J. There, it was noted that 
the interactions of a neutral fermionic field with the 
electroweak gauge-bosons of the SM are constrained by 
gauge-invariance to the magnetic dipole interaction 

£ M iDM = (^) X*^uB^x + c.c, (1) 
and the Rayleigh operators, 

£ RayDM = XX (cos x B^B^ + sin 9 x TrW^W") (2) 

R 

+ X15X (cos 9 x B^B^ + sin^TrW^W^) . 

R 

where in general A R =^ A R as we shall see below. Here x 
and x* are distinct Weyl fermions which may be two in- 
dependent Majorana fermions or part of a pseudo-Dirac 
pair. Ref. [2D] showed that in order to obtain the anni- 
hilation rate associated with the 130 GeV line reported 
in Ref. [IHl] the necessary dipole strength in the case 
of MiDM corresponds to a scale of about fi^ 1 ~ 2 TeV 
whereas the necessary Rayleigh scale in the case of 
RayDM is in the range A R ~ 500 - 600 Ge\Q Naive 
dimensional analysis then implies that the physics that 
resolves these non-renormalizable operators is somewhat 
strongly interacting (a point also made in [HI [20j [24]). 
In this letter we provide a concrete example of such UV 
physics with a simple renormalizable model that gives 
rise to a WIMP with both magnetic dipole and Rayleigh 
type interactions as featured in the MiDM and RayDM 
models. This UV completion helps to clarify the roles 
played by each of these interactions and their relative 
importance in different phases of the theory. 



Monochromatic gammas from the dipole operator in MiDM was 
first discussed by 28 . A more recent discussion where the early 
universe annihilation was suppressed by a large ~ 20 GeV split- 
ting (and a somewhat larger dipole) was recently considered in 



II. THE MODEL 

In addition to the WIMP state \ which is a Dirac 
fermion, we consider a messenger state, a Dirac fermion 
tp and a charged scalar tp, both of which are SU W (2) 
doublets with hypercharge Y = 1/2 and are heavier than 
the WIMP. They couple to the WIMP state through a 
Yukawa coupling which we denote by A. The Lagrangian 
for this model is given by 



X {%$ - mj x ~ 58m xCx + i>(0- Mf) ijj 



+ {D^tpY D^tp - M 2 s ^Lp + \tpx<p + h.c 



(3) 



where = 6 ^—igW^r" '—i\g' ' B ^ is the covariant deriva- 
tive associated with the SU W (2) x Uy(l) gauge-bosons, 
W° and respectively, and r a are the SU W (2) gener- 
ators obeying tr(r a r b ) = ^5 ah and related to the Pauli 
matrices through r a = \o a . Aside from its Dirac mass, 
m x , the WIMP states are split by a Majorana mass 5m. 

When the mass of the WIMP is much lower than that 
of the messengers, its interactions with light fields such as 
the photon and weak vector-bosons can be described by 
an effective Lagrangian. Gauge invariance forces these in- 
teractions to appear as dimension 5, magnetic dipole op- 
erator as well as dimension 7, Rayleigh operatorqj Since 
the model above is a renormalizable interacting theory 
these operators can be computed in perturbation the- 
ory. However, because we will be dealing with scenarios 
where the new states are not much heavier than the dark 
matter, it is important to include m /Mf corrections to 
these new operators (i.e., the form factors). In this let- 
ter we include all m /Mf effects at 1-loop order when 
computing the non-relativistic cross-sections relevant for 
phenomenology. 

We begin with the interactions of the WIMP with a 
single gauge-boson. These are generated through the di- 
agram shown in Fig. [T] Gauge-invariance forbids any 
coupling to the non-abelian SU W (2) fields and the most 
general vertex coupling to hypercharge consistent with 
Lorentz invariance can be written as, 



rV) = ^i(j J ) + <(f)^ 



1v F 2 (q 2 



(4) 



where the form- factors F±(q ) and i 7 ^*? 2 ) are given ex- 
plicitly in the appendijj^] The second part of this vertex 
corresponds to an effective dipole operator for the WIMP 
(if") X^ixvB^x with the dipole strength being 



I'x 



AY 
327T 2 Mt 



(5) 




FIG. 1. Magnetic dipole operator generated at 1-loop. 






(c) 



FIG. 2. The loop diagrams generating the RayDM operators 
at lowest order in perturbation theory. Diagrams (a), (b), and 
(c) represent two separate contributions where the external 
gauge-bosons are interchanged. 



where g' is the hypercharge coupling constant, q 2 is the 
momentum carried by the gauge-boson. More explicitly, 
the coefficient of the dipole operator is multiplied by the 
hypercharge and by the size of the SU W (2) representa- 
tion of the messengers in the loop, which in our case gives 
a factor of unity. Similar comments apply to the coeffi- 
cient of Fi(q 2 ). To lowest order in an expansion in the 
messenger mass these form-factors are 



Fi(q 2 ) = - 
F 2 (q 2 ) = 



y x q 2 hr 2 (3r 2 -3 - (2 + r 2 ) logr 2 ) ' 
6M f y (1-r 2 ) 2 
2r 2 (r 2 - 1 - logr 2 ) 
(1-r 2 ) 2 



2 After EWSB other, lower dimensional operators may appear in- 
volving the Higgs field, however those appear at higher loop order 
and are correspondingly much further suppressed. 

3 The F\(q 2 ) form-factor need not vanish as it is related to non- 
renormalizable terms of the form XJ 1 ' X^fv- Gauge-invariance 
only imposes the condition that Fi(q 2 ) should approach zero as 
q 2 -> 0. 



where r = Mf/M s . We include the effects of both Fx 
and Fi to all order in the messenger mass expansion in 
the cross-sections discussed below. 

The Rayleigh operators are generated by attaching 
another external gauge-boson to the loop diagrams, as 
shown in Fig. [2j In this case coupling to non-abelian 
gauge-bosons is possible as well. The Rayleigh scales as- 
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sociated with the non-abelian groups are 




( S?#Cf \ f 



J 7 , J 7 



(8) 



up to corrections of order O (m /M s ). (See for the 
equivalent calculation for a neutralino.) Here C/ is de- 
fined through the generators of the representation of the 
matter in the loop tr(i a i b ) = Cf5 ab and is equal to 1/2 for 
matter in the fundamental representation. In the abelian 
case this factor should be replaced by the square of the 
hypercharge times the size of the representation (so a 
factor of 2 x (1/2) 2 = 1/2 for the matter we consider). 
The functions jF{r) and J- are given explicitly in the ap- 
pendix. In general these are functions of all the masses 
in the problem as well as the Mandelstam variables s, t, 
and u. To lowest order in an expansion in the messenger 
mass they are given by 

r 2 (2 + 3r 2 - 6r 4 + r 6 + 12r 2 log(r)) 
Hr) = — 173 — > ( 9 ) 



(1-r 2 ) 3 
r 2 (3 - 3r 4 + 12r 2 log(r)) 
(1 - r 2 ) 3 



(10) 



with — and .F(l) = 1. When evaluating the 

resulting cross-sections we use the full contribution to 
1-loop order as discussed below. 

When a pseudo-Dirac state is present at freezeout the 
annihilation rate is dominated by the s-channel annihi- 
lation into SM charged fermion pairs through a j/Z ex- 
change. In the non-relativist limit appropriate for freeze- 
out this is given by, 



a v (XX -> ff) 



ira 2 q) (c w 1 F 1 + 2m x F 2 f 



m|/4m 2 



a / S w 2 + ( u / s w + ( 1 - m z/ 4m 



(11) 



Here c w (s w ) is the cosine (sine) of the Weinberg angle, 
qf is the electric charge of the fermion and Vf (a/) is the 
ratio of its vector (axial) coupling to the Z boson to its 
electric coupling. 

When only a single Majorana state is present the an- 
nihilation rate is dominated by the RayDM operators 
which contribute through WIMP annihilation into 77, 7Z, 
ZZ and W ± W IF [20]. By fixing the annihilation rate 
through these operators we can obtain a simple relation 
between the messenger mass scale and the Yukawa cou- 
pling A in the case of a WIMP with mass at 130 GeV, 



M mcd = 150 GeV (||) 



1/3 



av 

3 x 10- 26 ( 



1/6 
(12) 



The relation between the Yukawa coupling and the mes- 
senger mass in both scenarios is shown in Fig. [3j 
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FIG. 3. The Yukawa coupling vs. the messengers' mass fixing 
the freezeout rate at av = 3 x 10~ 26 cm 3 /s (solid) and av = 
6 x 10 _26 cm 3 /s (dashed). The upper (blue) curves are in the 
case when a single Majorana state is present at freezeout. The 
lower (black) curves are for the case when a pseudo-dirac state 
is present where the annihilation into SM fermions dominates. 

III. SCENARIOS FOR THE FERMI LINE 

With this in hand, we can explore the degree to which 
this UV scenario will actually provide the appropriate 
physics for the putative line signal observed by the Fermi 
satellite. We have essentially two scenarios to consider: 
MiDM, with a pseudo-Dirac state at freezeout and one 
with a pure Majorana state (where the heavier dark state 
is decoupled). 

To begin with, it is worth emphasizing that except for 
the most strongly coupled region of parameter space for 
the MiDM scenario, the 77 signal in the present day is 
dominated by the dimension-7 Rayleigh operators, as ar- 
gued in |20) . where the hybrid scenario was proposed. 
Consequently, for both MiDM and RayDM scenarios, 
the ratios of 77 and jZ are determined essentially by 
the representation of the messengers. For doublet mes- 
sengers, the re lative si zes of the Rayleigh operators is 
cos# x = gy/y 9y + 9 ~ 0.29 (i.e., the sizes are deter- 
mined by the gauge couplings alone). 

For a doublet, this yields a ratio of ct 77 / (cr 7 ^/2) ~ 2.2 
and <7t i/(cr 77 + \/2a 1 z) — 5. Thus, it is a natural ex- 
pectation for models dominated by Rayleigh annihilation 
from a doublet loop that both lines should be visible. 
Secondly, we emphasize that there is no issue with con- 
straints from continuum emission in these models. The 
ratio of monochromatic to continuum photons is safely 
below the limits of [TT1[13] . This is true irrespective 
of whether the freezeout occurs through the dipole or 
through the Rayleigh operator (as the continuum ratio 
is the same in both cases). Essentially, continuum emis- 
sion is a problem for models where the line annihilation 
is produced at a higher order in some perturbative ex- 
pansion from the continuum. Here, they are produced at 
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FIG. 4. The solid (dashed) curve depicts the annihilation rate 
XX — * 77 ( — * 7^) as a function of the mass of the messenger 
in the case of a pseudo-Dirac WIMP. The Yukawa coupling 
is fixed by requiring thermal freezeout with an annihilation 
cross-section of av = 6 x 10 _26 cm 3 /s. 



the same order and the lack of any a priori problem is 
obvious. 

Let us focus for a moment at the MiDM scenario. Here, 
the freezeout occurs via the dipole annihilation into //, 
while the present day annihilation is dominated by the 
Rayleigh operator. For messengers in the 150-200 GeV 
range, the coupling a « 1 is strong, but still perturba- 
tive, and the theory is calculable. The annihilation into 
photons (crv) 7J + ^(crv) jZ rj 10 _27 cm 3 s _1 is precisely the 
right value to explain the observed excess. Intriguingly, 
the dipole here is slightly smaller than 10 -3 fijy, which 
is the size previously argued to explain the DAMA an- 
nual modulation result [30 . Thus, it is conceivable (and 
is already strongly constrained from direct detection ex- 
periments) that such a scenario could also yield an ex- 
planation of the DAMA result. 

In the case where we have only a single Majorana 
fermion (both in the present universe as well as at freeze- 
out), we must have a truly strongly coupled theory to 
generate the Rayleigh operator of the appropriate size. 
Again, we have both B^B^ and W^W^ operators 
with cos9 x — 0.29. Assuming that together, these pro- 
vide the appropriate relic abundance we have a Rayleigh 
scale of ~ 550 GeV (a difference of 2 1 / 6 from the Dirac 
case). Here, normalizing to freezeout, we expect a cross 
section of er 77 -I- l/2a 7 z ~ 6 x 10~ 2 'cm 3 s _1 , which is 
large, but perfectly acceptable for a slightly flatter halo. 



IV. CONCLUSIONS 

The recent evidence for a possible dark matter sig- 
nal in gamma rays has prompted a reexamination of 
the interactions of dark matter with light. For a Ma- 



jorana fermion, the leading operators are a dimension-5 
dipole operator in the presence of a nearby excited state 
(the MiDM scenario) or a dimension-7 Rayleigh opera- 
tor more generally. The scales of the operators (~ TeV 
for the dipole and ~ 600 GeV for the Rayleigh operator) 
suggest that the UV completion is at or near the weak 
scale. 

In the presence of the simplest possible theory that 
generates these - namely, a loop of electroweakly charged 
messengers - we can understand the overall phenomenol- 
ogy for the 130 GeV line and more generally. We have 
found that for most of the parameter space, except- 
ing only the most strongly coupled, the Rayleigh op- 
erator dominates the present-day 77 signals, while the 
dipole annihilation \X* ~~ ^ ff (which is inaccessible to- 
day) dominates freezeout. With weak-scale messenger 
masses and a strong, but perturbative (a ~ 1) coupling, 
MiDM provides a natural framework to explain these sig- 
nals, without any apparent conflict from the data. If 
the charged matter carries transforms under some strong 
gauge group such as in e.g., Sister Higgs models [3T], such 
couplings are not unreasonable. Intriguingly, the gener- 
ated dipole is also of the size necessary to explain the 
DAMA modulation. Without the excited WIMP state 
for freezeout, annihilation through the Rayleigh opera- 
tor also provides a viable scenario both for the 130 GeV 
line and relic abundance, but at the cost of both rel- 
atively light matter and very strong (non-perturbative) 
couplings. 

In both scenarios the strength of <7 77 /(|<7 7 2) is de- 
termined entirely (at leading order in weak couplings) 
by the 517(2) x U(l) representations of the matter in 
the loop. For (2, ±1/2) messengers, the ratio is roughly 
2.2:1, consistent with recent suggestions. 

The low scales of the new matter imply that the effec- 
tive operator approach is limited in its quantitative ap- 
plications. Indeed, including all orders in the m x /M mess 
expansion tend to enhance the annihilation rates both 
in the late and early universe. Nonetheless, normalizing 
to the appropriate relic abundance, the present day sig- 
nals are not dramatically changed when including these 
effects, only their interpretation in terms of the masses 
and couplings of the new states. 

In summary, one can understand the effective theory 
of MiDM and RayDM with a simple UV completion that 
gives the relative signals and scales in different regions of 
parameter space. Remarkably, this simple completion is 
better than just a toy model, providing a successful the- 
ory at perturbative coupling. One must accept a some- 
what strongly interacting theory, but there is no reason 
to think our calculational challenges prohibit nature from 
realizing such scenarios. Should the Fermi line prove to 
be true evidence of dark matter, this model may help 
provide qualitative and quantitative understanding of the 
signal. 
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Appendix A: Calculating the Dipole and Rayleigh Coefficients and Form Factors 



In this supplement we provide explicit and detailed expression for the results quoted and used in the paper. We 
also offer simplified expression in some limits of physical interest. 



1. Single gauge-boson vertex 



The amplitude with the gauge-boson attached to the scalar is 



9,M 



V \AA/ 



Pi. 



v{jP2) («A) 



d±k + Mf) I 

(2tt) 4 (pi - k) 2 - M) k 2 - Mf 



(-ig(2k - qT) 



(k - q y - mi 



(-i\)u(p 1 ) 



2X 2 v(p 2 ) fcPx^-f ^ 2 ^ ' 



2 ' (47r) d / 2 A^ d/2 

i 2m x z ((1 - z)m x + M f ) (Y + ia^q v /2m x ) 
+ 32^ /T 



«(pi). 



(Al) 



Here d 3 x = dx dy dz are the usual Feynman parameters, restricted toi + j/+z = l and we used dimensional 
regularization with e = 4 — d to define the integral. The denominator is given by 



A s = zMf + (x + y)M 2 - z(l ~ z)m 2 x - xyq 2 



(A2) 



The divergent piece in Eq. (Al) cancels against the divergent piece in the second diagram where the gauge-boson is 



attached to the fermion line. The amplitude for this process is 

iM f - 




d 4 k 



i {-f -£ + M f ) i(?i-t + M f ) 

A (isrf) — 



(2tt) 4 (-pa - k) 2 - M 2 v " a ' ' (pi - k) 2 - M 2 k 2 - M 2 

r(2-|) i 



(-iX)u(p 1 ) 



2x 2 v( P2 ) / d^ 2 r (i - f) 



(Air) d / 2 A 2 f r d/2 
i {{zm x + Mf) 2 + xy q 2 )Y - 2i x (zm x + Mf) a^ u q v 



32tt 2 



)<Pi). (A3) 



with 



Af = zM 2 + (x + y)M 2 - z(l - z)m, 2 - xy q 2 



(A4) 



() 



combining both diagrams the limit e — > yields a finite expression. The result can be written in terms of two 
form-factors as in Eq. (3) in the text 

FV) = rf Fl ( g 2) + i (bC) <J^q„F 2 (q 2 ) (A5) 

The dipole strength fj, x is given by 

g'A 2 

327r 2 M/ : 

and the form-factors are 

F,(,r:,„ N \I f .MJ = ^ I ^ \ l + - ^— ■ ' A J ■ ' ' lo,(^]) ( A7) 



and 



a'A 2 f o f 2m (Mf + (l-z)m ) (M f + zm )"+oto 2 /A,\ 

W;".„m„m,) = ■?„ a ^ J - ^ ^ ^ + bg ( 

When the momentum exchange is small q 2 — > the charge form-factor approaches zero as it should since the WIMP 
state is uncharged 

F l{ q 2 ) ^ _/^ 2 ( 3, 2 -3-(2 + , 2 )log, 2 ) , v 

v ; 6M/ (l - r 2 ) \ f J y ' 

where r = Mf /M s . Similarly, the dipole form-factor can be simplified in the limit of heavy messengers 

. , s Mf,M»»g 2 ,m v 2r 2 (r 2 — 1 - log r 2 ) f ,\ 

F 2 (q 2 ) -L K — + O (My 1 ) . (A10) 

(1 — H) v y 

More importantly are the simplification that occur in the non-relativistic limit relevant in the case of the annihilation 
into SM fcrmion pairs, XX ~^ f f where q 2 = 4m 2 to lowest order in velocity. For simplicity we only consider the case 
of equal messenger masses since the most general case yields extremely complex and unilluminating formulas. The 
charge form-factor is 



g'X 2 



« 2=4m x 607r 2 m|,/l -m 2 



( 2m\ (1 - m 2 ) 3/2 + (-16 + 18m 2 - 3m* + m\) arcsin (m x ) (All) 
+ (4 - m 2 ) 2 ^4 - 5m 2 + m 4 , arctan TO * J ) 



4 — m 2 



« 2=4m x 15^6^4-???, 



X 



( 2^Jl - m 2 ^4 - m 2 (-16 + m 2 (12 + m x (15 + 4m x ))) arcsin (m x ) (A12) 



+ {m 3 x yjA -m 2 (4 + m 2 ) + 2(-2 + m x )(2 + m x ) 2 (-8 + m x (4 + m x (7 + 4m x ))) arctan ( 



4 — m 2 



where fh x is the mass of the WIMP in units of the messenger mass. Since the WIMP is light than the messengers a 
power expansion in m x is likely more useful and illuminating 



F, 



( <A 2 \ { _ 2 19mt 19m x 29m» \ _ oN 

^ = - (M r* + v + v + v j + ° k- io ) (Ai3) 

= l + ^3 + ^ + 3 ^+0(m-) (A14) 

g2=4m2 3 2 45 10 Vx; V7 

In Fig. [5] we plot the WIMP's dipole strength generated from the loop as a function of the coupling to the heavier 
charged states. Evidently, the coupling needs to be sizable in order to obtain a sufficiently large dipole even for 
charged states which are only slightly heavier than the WIMP. 
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FIG. 5. The WIMP's dipole strength in units of the nuclear magneton fi N as a function of its coupling to the heavier charged 
states. For a given choice of the coupling, the messenger mass is set by requiring the annihilation rate to SM fermion - anti- 
fermion pair to be equal to 3 x 10 _26 cm 3 /s. The curve asymptotes at large coupling since at that point the annihilation rate 
is dominated by the dipole contribution. 



2. Two gauge-boson vertex 

There are a total of seven loop diagrams contributing to the two gauge-boson vertex, three of which are related 
by a simple exchange of the two gauge-bosons. The first diagram where both external bosons are connected to the 
fermion line is the most complicated one and is given by 



Pi, n pi, v 




where 

k = k! + (xip! - X2P2 + x 3 (pi - p 3 )) (A16) 
Ai = Alj(xi + X2 + X3) + MgXi — SX1X2 — tx 3 — mfxix 3 + tx\Xy, + m\x2X 3 (A17) 
- sx 2 x 3 - ux 2 x 3 + tx 3 + m 2 x {x\ — £Ci(l - 2x 2 - a; 3 ) - x 2 (l - x 2 - 3x 3 )) 

Here s, t, u are the usual Mandelstam variables (s = (pi + P2) 2 , t = (pi — p 3 ) 2 , and u — (pi — P4) 2 ) and m 3 4 = p 3 4 
are the external bosons' mass, which we keep explicit as a check on the calculation since they should drop out when 
extracting the final answer for the coefficient of the Rayleigh terms. Shifting the momentum to fc', the momentum 
integral is finite and easily doable. The second amplitude, iA^2, has the same form but with the two external bosons 
exchanged, 



iM.2 = O v, m 3 <-} rri4,p 3 ■<-» p^, t o u) 



(A18) 



The next diagram has one external gauge-boson attached to the fermion line while the other to the scalar line 

Pi," 



^3 = /' 




•*r.\-.-(M / < / ( fc , 2 _ A3) 4 l«(Pi) (A19) 



where, 



k = k' + (xipi +x 2 (p3 -pa) + ^3) (A20) 
A3 = Mj(xi + X2) + M s (X3 + X4) — UX2 + TO3X1X2 — SX1X2 — tx\X2 + ux\ — m 3 x 3 (A21) 
+ 7713X1X3 — txix$ + 7773X2X3 + UX2X3 + to 3 xJ + m 2 ^ [x\ — x 2 a;3 — X\{1 — 3x2 ~ ^3)) 

The fourth amplitude, zA^4, is similar but with the two external gauge-bosons exchanged. The next diagram is 

Pi - k 

Pi Pi 

iM 5 = — » > » — ] — * — 

k //Pi +P2-k 

( f . f d*k (? 1 -fi + M f )(2k-p 3 r(2k-p 1 -p 2 -p 3 V\ 
= -3!^)(/^/^ ^ (fc , 2 _ As)4 )«<») (A22) 

where, 

k = k' + (x x pi + x 2 (pi +p 2 ) + ^3) (A23) 
A5 = MfXi + M 8 (x2 + X3 + X4) — SX2 + SX1X2 + sx 2 — 7TI3X3 + 777,3X1X3 — txiX3 (A24) 

+ 2777,3X2X3 — <X2X3 — UX2X3 + 7773X3 + 777^ (x\ — Xl(l — X3) + 2X2X3) 

The sixth amplitude, iAAq, is similar but with the two external gauge-bosons exchanged. The final diagram is 

Pi - k 

Pi > < P'i 
1M7 = * * * 1 * 



feV^ ^jPi+n-k 



Pi,V 



where, 



= 2 2! " 2 < V2 "^( / / (Hi ' ( fc >2-A 7 )3 j U(Pl) (A25) 



fc = /c' + (xirJi+x 2 (pi+P2)) (A26) 
A7 = MfXx + (X2 + X3) — m (1 — xi)xi — SX2 + SX1X2 + sx 2 (A27) 
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FIG. 6. On the left pane is a plot of the relative strength of the coefficients of the two Rayleigh operators as a function of the 
ratio of the fermion to scalar messenger. On the right pane is a plot of the F£ form-factor as a function of the messenger mass, 
assuming m x = 130 GeV and messengers of equal mass. Fg asymptotes to unity in the large messenger limit. 



Expanding the amplitudes in powers of inverse messenger mass the sum of the diagrams can be shown to vanish 
up to order O (M~^ ss ) . At this order the first contribution to the Rayleigh operators appear. The matching between 
the coefficient of the terms that appear in the Lagrangian and the terms in the amplitude is as follows, 



XXFuuF^ 



and 



=> Av(p 2 )u(p 1 ) (p 3 -Pig^ - P3P l {js fl (p 3 )s l ,(p4) 

(P2)( - ^p CTK A7 P 7 ff 7 K 7 A )«bi)(4^ Ql/ ^3^)e M (p 3 )£,b4) 



Here £^{pz) and £ u (pa) are the external gauge-boson's polarization vectors. In the large messenger mass expansion 
the Rayleigh scales are found to be 



1 



A 3 , 



g 2 A 2 
48M|tt 2 



2\2 



9^ 



\ 48M]it 2 



with 



Hr) 
Hr) 



r 2 (2 + 3r 2 - 6r 4 + r 6 + 12r 2 log(r)) 

(1 -r 2 ) 3 
r 2 (3 - 3r 4 + 12r 2 log(r)) 
(1 -r 2 ) 3 



(A28) 
(A29) 

(A30) 
(A31) 



Here r = Mf /M s is the ratio of the scalar messenger's mass to that of the fermion messenger. The form-factors obey 
J-{1) — and .F(l) = 1. A plot of the relative strength of these two coefficients is shown in Fig. ^ below. 

In the case of non-relativistic annihilation these results are inadequate when the messenger mass is not much larger 
than the WIMP mass because the kinematical variable s = Am 2 is of the same order (or larger) than the messenger 
mass. Not only is the expansion in inverse powers of the messengers mass is inadequate, one must worry about higher 
dimensional operators that could contribute. However, as we now argue, to one-loop order in perturbation theory 
such higher order operators do not contribute. In other words, no other terms in the amplitude aside from those 
associated with ixi5XF^F fil ' contribute to the non-relativistic cross-section. 

We first note that any terms with pi j2 ■ e(P3,4) are velocity suppressed and can be neglected since the incoming 
momenta (which are mostly time-like) are contracted with the outgoing polarizations (which are space-like). Any 
terms with p 3 ■ e{p 3 ) and pn ■ s{pi) vanish identically because of the transversality of the gauge-boson's polarizations. 
Therefore, by momentum conservation, p 4 • e{pz) and p% • e{pi) are both velocity suppressed and can be neglected. 
Thus, any term with external momenta contracted against the gauge-bosons' polarization can be neglected. 
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Aside from the terms associated with the dual Raleigh operator, ixisX^iyF^ , the only surviving terms are ones 
associated with ixxF^F^ which is itself velocity suppressed. So we need only extract the coefficient of the dual 
Raleigh operator albeit to all orders in the messengers mass. For simplicity we quote the result in the case of equal 
messenger mass 



A 3 \ 48M*n 2 

R \ 7 , 



F 3 (s,t,u) (A32) 



JL__ / g 2 X 2 
where the form-factor F 3 is given by the integral 

f 4a 

rh x (l - rh 2 , (— x 4 + (2xi + X4)(2x 2 + X4))) 



, , 2 ^ / ,4 , 2(2rh x + 4x 3 +m 2 x (2x 4 + (l-x 3 )(-3x i + 2(2x 1 + x 4 ){2x 2 + x i )))) 
F 3 ( 4 "\,- TO v ,- m v ) = I d x [ — j± — — — (A33) 

(A34) 



m x (l + m 2 (xi - x 2 )(x 3 - x A )) J 
The enhancement due to this form-factor is plotted on the right pane of Fig. [6j 

3. Cross-Sections 

The non-relativistic annihilation cross-sections associated with the two gauge-boson vertex (RayDM) were taken 
from ref. [20] and are reproduced here for completion. The general expression is given by 

a 2 m 4 

<T(XX^VV)V = ^ 7 +1C VV (A35) 

R 

with the kinematic functions fC vv and couplings g vv defined as 

/C 77 = 1 g 7 j — cos 6 x cos 2 &w + sin 9 x sin 2 9w (A36) 

/C 7 z = 2 M - ff7Z = \ (sm6 x - cosflj sin(20 w ) (A37) 

( 777 ^ \ 

JC ZZ = 1 - -f I g zz = cos 9 x sin 2 9 W + sin 9 x cos 2 W (A38) 



3/2 



K W w =2 I I g ww =sm9 x (A39) 



In the case of messengers in the (2, |) representation of SU W (2) the angle cos# x = g' 2 / \J g' 4 + g 4 = 0.29. 

The most important annihilation cross-section associated with the single gauge-boson vertex is that of the WIMPs 
into SM fermion pairs. If present, this annihilation mode dominates over all other channels. Using the general 



expression for the one gauge-boson vertex, Eq. ( A5 1 we obtain the following matrix element 

\ E W - 2 / 4g ' 2 , 2 2 ( Fl (2m 4 + .s 2 - AmH + 2st + 2t 2 ) + F X F 2 U^s 2 (*£) c w ) (A40) 

+ Fl (-s(ml - 2m 2 (s + t) + t(s + t)) (^) * + S (-m 4 + 2m 2 (s + 1) - t(s + t)) (^) * (c w 2 - s w 2 )) 

x (mf - 2m\s + s 2 (l + a) + vj) + (m z - 2m 2 z s + s 2 (l - a) - vj)) (c w 2 - s w 2 ) + 4s (s - ml) w/s w c w ) 
where s and t are the usual Mandelstam variables s = (p\ +P2) 2 and t — (pi — p 3 ) 2 ■ 
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